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Read in B-Rep model M 
Assign Face IDs F to all faces in M 
Assign EdgelDs E to all faces in M 
Set voxel size S 

Find B-Rep model bounding box 

Find discrete bounding box coordinates X, Y ,Z 

Create voxel text file T 

for i in X in increment of S: 
for j in Y in increment of S: 
for k in 2 in increment of S: 



construct voxel v of size S using AC IS 
test intersection of v with M using ACIS 
if v intersects with M: 
store "1" in T 

find faces intersecting with v using ACIS 
find edges intersecting with v using ACIS 
store F and E for all intersecting faces in v 

else: 

store "0" in T 
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A step of converting tubular structure to prismatic structure 

1 . Voxelize a part volume 

2. Voxelize a cavity volume 

3. Gain internal cavity voxels using a degree of visibility 

4. Merge part voxels and internal cavity voxels 
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FIG. 35: (a) 3D model, (b) Levels of detail (LOD) for n = 2, (c) LOD for n = 3, (d) LOD 

forn = 4,(d)LODforn = 5 
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• start() // Find die Starting voxel 
{ 

- Find a terminal node 

- Else choose top:left:back most voxel 

} 

• marchQ 

- for each next voxel 
" { 

• findNeighbors() 

• If(nNeigh = Nil) // nNeigh - Number of Neighbors 

- NODE; Terminate 

• If(nNeigh — NORMAL Neighbors) 

- Change voxel content to currentEntitylD 

- Push voxel into current entity's stack 

- Iterate 

• If(nNeigh > NORMAL Neighbors) { 

- Check isNode() for all neighbors 

- Find the node and step into it 

- Push the current node to the Node stack 

- For all branch { 

march along() 

" } 

• findLoopEntities(){ // Simple Loops 

Trace back the previous nodes 

Add the corresponding entities to the loop 

• Pop node from the current node stack 

• iterate 
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(ii) Results after refining the query in (i) using moment invariants (Second Step) 

FIG. 63 



